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ABSTRACT 

We employ Schwarzschild's method of orbit modeling to constrain the mass profiles 
of the central lens galaxies in Q0957+561 and PG 1115+080. We combine the measured 
central projected stellar velocity dispersions of these galaxies with the self-similar radial 
profiles of the rms velocity and of the Gauss-Hermite moment h.4 observed in nearby 
galaxies for <^ R <^ 2i? e ff- For Q0957+561, we find a 16% uncertainty in the galaxy 
mass, and formal 2-a limits on the Hubble constant of Hq = 61 IjJ km s _1 Mpc -1 . For 
PG 1115+080, we find that none of the viable lens models can be ruled out, so that Hq 
is not yet strongly constrained by this system. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: structure — galax- 
ies: fundamental parameters — galaxies: kinematics and dynamics — gravitational 
lensing — quasars: individual (Q0957+561) — quasars: individual (PG 1115+080) - 
distance scale — dark matter 

1. INTRODUCTION 

The gravitational lens system Q0957+561 (Walsh, Carswell, & Weymann 1979) has been 
modeled extensively in an effort to determine the Hubble constant Hq from measurements of the 
time delay between its two primary images. Since the long-running dispute over the time delay 
measurement has been resolved in favor of the short delay (Schild & Thomson 1997; Haarsma et al. 
1997, 1998; Kundic et al. 1997), the largest remaining uncertainty arises from the mass model. The 
lens consists of a cluster with a large central elliptical galaxy Gl. The asymmetric radial positions 
of the images accurately constrain the parameterized radial mass distribution of Gl (Grogin & 
Narayan 1996, hereafter GN). However, the presence of the cluster introduces a degeneracy in 
the overall mass normalization of Gl (Falco, Gorenstein, & Shapiro 1985; Gorenstein, Falco, & 
Shapiro 1988), and thus in the determination of the Hubble constant (#0 <x °~o> where a$ is a 
velocity dispersion characterizing the mass of Gl). Therefore, additional independent constraints 
are needed on the relative contributions of Gl and the cluster to the 6" image separation. These can 
be obtained by inferring the mass of the cluster from cluster dynamics (Garrett, Walsh, & Carswell 
1992; Angonin-Willaime, Soucail, & Vanderriest 1994), from hot intracluster gas X-ray emission 
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(Chartas et al. 1995, 1998), or from the weak lensing of background galaxies (Dahle, Maddox, & 
Lilje 1994; Fischer et al. 1997). Alternatively, one can infer the mass of Gl from stellar dynamical 
measurements. 

The quadruple lens PG 1115+080 (Weymann et al. 1980) is the second system with a well- 
determined time delay (Schechter et al. 1997; Barkana 1997b). The projected mass of the primary 
lens galaxy G inside the ring of images can be precisely determined (unlike the case of Q0957+561, 
uncertainties in the mass distribution of nearby galaxies have only a minor effect on the models). 
However, the geometry of the system does not permit distinguishing between different mass profiles. 
This has important consequences for the Hubble constant, since Ho can vary by 40%, depending 
on the mass model assumed (see Keeton & Kochanek 1997; Courbin et al. 1997; Saha & Williams 
1997; Impey et al. 1998). Stellar dynamical measurements of G may be helpful for breaking the 
degeneracy in the mass model, and thus in Hq. 

Falco et al. (1997) measured the central projected stellar velocity dispersion of Q0957+561 
Gl to be <7p = 279 ± 12 km s _1 , improving on an earlier measurement by Rhee (1991). Similarly, 
Tonry (1998) measured the central dispersion of PG 1115+080 G to be a p = 281 ± 25 km s _1 . 
However, the conversion of the measured a p to oq is subject to systematic uncertainties, which 
include the unknown anisotropy structure of the stellar orbits, the radial variation of the mass-to- 
light ratio, and the ellipticity of the galaxy. Previous galaxy models have arrived at a relatively 
small uncertainty in this conversion by making arbitrary simplifying assumptions. Kochanek (1993, 
1994) considered a singular isothermal mass model with constant anisotropy. For Q0957+561, GN 
used these dynamical models, but limited them to be nearly isotropic, leading to a 2% systematic 
uncertainty in Oq (and thus in Hq). Barkana (1997a) also assumed near-isotropic orbits, leading to 
an uncertainty in Ho of 4%, but he noted that allowing for more anisotropy gives an uncertainty of 
14%. In fact, dynamical studies of nearby elliptical galaxies have clearly demonstrated that there 
is little basis for the assumption of isotropy, or even of constant anisotropy. 

Binney & Mamon (1982) and Tonry (1983) first illustrated that an elliptical galaxy's surface 
brightness and projected stellar velocity dispersion profiles, I(R) and o~ p (R), could not determine 
both its mass distribution and its anisotropy profile. Richstone & Tremaine (1984) and Katz 
& Richstone (1985) used orbit modeling methods to demonstrate that the conversion from the 
dispersion profile cr p (R) to a mass parameter <Tq can be uncertain by an order of magnitude. Further 
theoretical studies (Dejonghe 1987; Merritt 1987; Merrifield & Kent 1990; Gerhard 1991; Dejonghe 
& Merritt 1992; Merritt k, Saha 1993; Merritt 1993) demonstrated that complete knowledge of 
the stellar line-of-sight velocity distribution (LOSVD) C(v p ,R) gives a unique solution for the 
two-integral distribution function (DF) f(E,L) in a known spherical potential <E>, and may even 
strongly constrain an unknown $. The general efficacy of incomplete knowledge of C is less clear, 
but important constraints on / and <E> could be further provided by large-radius measurements of 
higher-order velocity moments (e.g., the Gauss-Hermite moments hi — van der Marel & Franx 
1993; Gerhard 1993; see also Rix & White 1992; Zhao & Prada 1996). Rix et al. (1997) and 
Gerhard et al. (1998) made dynamical fits to nearby galaxies with higher-order moments (h^jh^) 
measured to ~ 2.5i? c g; they determined the total mass to ~ 10-15%, and ruled out a constant 
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mass-to-light ratio with > 99% confidence. For the galaxy NGC 2434, the solutions typically had 
nearly constant radial anisotropy (3(r) = (1 — Vg/v^) ~ 0.5, while NGC 6703 showed an anisotropy 
rising from (3 ~ 0.1 at the center to ~ 0.4 near i? e flf- Thus the assumption of f3(r) = 0, or even of 
constant (3{r), is certainly unwarranted and probably incorrect. 

While our current knowledge of elliptical galaxies does not permit us to make arbitrary as- 
sumptions about the anisotropy of the DF, we can impose the constraint that the unmeasured 
cr p (R) and h^(R) profiles of these lens galaxies are similar to those of other galaxies (provided 
they are universally homologous). In this study we rigorously consider the utility of the central 
<7p measurement for determining Hq from these two lens systems. To ensure physically correct, 
robust results, we use a spherical orbit modeling method after Schwarzschild (1979), Richstone 
& Tremaine (1984), and Rix et al. (1997) — a fully general way to construct realistic models of 
a galaxy, given an assumed potential. In §2 we review the observational constraints on the lens 
galaxies, and introduce constraints on their LOSVD profiles by demonstrating that the profiles of 
nearby elliptical galaxies are self-similar. We describe our modeling method and demonstrate it 
with a test-case problem in §3. In §4, we report the range of model solutions for Q0957+561, and 
discuss the implications for Hq. We examine the solutions for PG 1115+080 in §5. In §6 we give 
our conclusions. 



We first review the observational data for Q0957+561 Gl and PG 1115+080 G in §2.1. As 
discussed in §1, a simple measurement of the central projected stellar velocity dispersion a p cannot 
strongly constrain the mass of a lens galaxy, necessitating further a priori constraints. However, 
rather than making unjustified assumptions about the galaxy's anisotropy, we impose conditions 
on its observable properties by requiring it to have an LOSVD profile consistent with those of 
better-observed galaxies. To this end, we extract "mean profiles" of v Tms (R) and hi(R), including 
their galaxy-to-galaxy scatter, from the observational data available for nearby elliptical galaxies^]. 
If there proves to be a universal shape to these profiles (i.e., there is little scatter between galaxies), 
we can use them as additional constraints (see §2.2). 



As introduced by van der Marel & Franx (1993) and Gerhard (1993), C is parameterized by the Gauss-Hermite 
moments, 



where w = (v p — v p )/a p , 70 is the line strength, (7, v p ,a p ) are the coefficients for the best Gaussian fit to C, and 
Hi(w) are the Hermite polynomials. A perfectly Gaussian profile will have hi = 801 and a p — a p . For spherical 
systems, the fourth moment /14 is a useful constraint on the DF; its value is typically positive (peaked LOSVD) if 
/ is radially anisotropic, and negative (flat-topped LOSVD) if tangentially anisotropic. Strictly speaking, since we 
apply the data to non-rotating models, we should use the Gauss-Hermite moment Z4 corresponding to the even part 
of the LOSVD (van der Marel et al. 1994 §5.1), but such data are not available, and the correction is probably small. 
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OBSERVATIONAL CONSTRAINTS 
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2.1. Observations of Q0957+561 and PG 1115+080 

Accurate observations of Q0957+561 Gl are difficult because of the nearby bright quasar 
image, but measurements have been made of the stellar surface brightness profile I(R) and the 
central velocity dispersion a p . Bernstein, Tyson, & Kochanek (1993) measured 1(a) along the 
major axis (a = 2'.'5-18"3) in the i?-band, and Bernstein et al. (1997) measured it at smaller radii 
(a = 0"l-5"7) in the F-band. The galaxy is rather round, with an ellipticity increasing with radius 
(e ~ 0.1, 0.2, 0.4 for a ~ O'.'l, T'5, 10"). We combine the R and V data, assuming a simple offset of 
1.34 magnitudes, and map the data to the intermediate radius m = \J a(l — e), resulting in a profile 
I(m) with m = 0'.'1-12"8, where R e s ~ 4'.'5. To ensure a reasonable I(m) profile, we also impose 
weak constraints at large and small radii, resulting in 50 data points over the range m = 0'.'0-29"4. 
Falco et al. (1997) measured the central dispersion of Gl to be a p = 279 ± 12 km s _1 inside a 0'.'6 
(1.9/i — 1 kpc) radius. Their data suggest a rise in the dispersion at the center, with a p = 316 ± 14 
km s _1 inside 0'.'2, and a p = 266 ± 12 km s _1 outside 0'.'2. Since such a rise is inconsistent with 
the radial dispersion profiles of nearby galaxies and with the subsequent dispersion measurements 
of Tonry & Franx (1998), we adopt the total binned measurement of 279 km s -1 . 

Grogin & Narayan (1996) modeled the lensing properties of Q0957+561, using two different 
parametric mass models for Gl. We will compare our results with their softened power-law sphere 
(SPLS) model results, where their density profile is p(r) = po(l + r 2 /r 2 )~ a / 2 . Their best-fit 
parameters are (a = 1.92±g;g| J r c = 0'.'058l^^|, a E = 2^40±g;i§|, at 2 a), where the deflection 
parameter «e is related to the central density po- This best-fit model is a poor fit to the data 
(x 2 per degree of freedom of 6.9), and uses a position for Gl that has been shown to be incorrect 
(Bernstein et al. 1997). But for want of a reanalysis of the lensing constraints, we will adopt this 
solution in our models. 

Impey et al. (1998) found that the galaxy PG 1115+080 G is nearly circular, and fit well by 
a de Vaucouleurs profile with i? e fj = 0'.'59 ± 0"06. We approximate the surface brightness of this 
galaxy by a Hernquist (1990) model with a break radius of a = 0'.'325 (i? e flf — 0"59), modeled as a 
profile I(m) with 21 bins over the range m = 0'.'0-2"l. Tonry (1998) measured the central dispersion 
of G to be <Tp = 281 ± 25 km s" 1 inside a 0"6 (1.7/i _1 kpc) radius. Impey et al. (1998) fit the 
lensing constraints with three standard mass models for G: a singular isothermal sphere (a = 2, 
r c = 0"), a modified Hubble profile (a = 3, r c = 0'.'2), and a constant mass-to-light ratio model. For 
all the models, the total projected mass inside R = 1"15 was found to be (1.24-1.39) xl0 11 /i _1 M Q , 
depending on the mass model assumed for a nearby galaxy group. 



2.2. Self-Similarity of Kinematic Profiles 

Bender, Saglia, & Gerhard (1994) measured the LOSVD out to ~ R c d for a large, "unbiased" 
sample of galaxies, and derived the profiles v p (R), a p (R), hs(R), and h^R). For our data set, we 
take from this sample 80 profiles from 28 elliptical galaxies. We find the rms projected velocity -u rms 
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by numerically integrating the positive line profile C(v p ) with fit parameters {v p , a p , h±} (see 
van der Marel & Franx 1993 §2.4). To make a scale- free comparison of the profiles, we renormalize 
to the central v rms (inside 0.14i? c g, for direct comparison with the central a p measurement of 
Q0957+561 Gl; or inside l.lR eS for PG 1115+080 G), and rescale radially by R cS . We use both 
major and minor axis data, mapping them to the intermediate radius m. We combine the data 
in radial bins, spaced such that the number of points in each bin is nearly constant (~ 36). To 
improve the large-radius constraints, we divide the last bin (0.7-1.8 R c r ) into three bins, with 12-13 
points in each bin. As shown in Figure lb for the case normalized to Q0957+561 Gl, the resulting 
(rescaled) "mean v rms profile" is nearly constant (t> rms ~ 0.9 ± 0.1 at 1.5i? e ff), while the dispersion 
a p = ((v p — Vp) 2 ) 1 / 2 decreases with radius (a p ~ 0.6 ± 0.1 at 1.5i? c ff ) due to increasing rotational 
support in the outer parts of the low-luminosity galaxies (see Davies et al. 1983; Fisher, Illingworth, 
& Franx 1995). We thus find a remarkably universal, flat rms velocity profile inside 1.5i? e ff, for 
all elliptical galaxies regardless of other properties. Note that v rms is not only less variable than 
<7p, but is also a better physical probe of the gravitational potential. The corrections in the mean 
profiles due to h% and h± are small, so that a v ~ a p to ~ 0.3% accuracy, and v rms ~ (v p + cr p ) 1 / 2 
to ~ 5%. 

To produce the mean h± profile, we use the same binning procedure as for v rms (see Figure 
la). We also add data from a sample of six galaxies (Carollo et al. 1995; Statler, Smecker-Hane, 
& Cecil 1996; Gerhard et al. 1998) for which /14 has been measured at larger radii (1.7-4.9i? e ff); 
we combine these points into three additional bins, with ~ 12 points in each bin. Note that these 
points are not a large sample, and comparisons of different authors' data sets suggest that the 
errors in large radius velocity data are generally underestimated. The mean and dispersion range 
from h 4 ~ 0.002 ± 0.03 in the central bin to 0.02 ± 0.06 at l.5R cS , and the LOSVD is everywhere 
consistent with Gaussianity (/14 = 0). 

Because our mean profiles are derived from a general data set of elliptical galaxies, they may not 
accurately represent a brightest cluster galaxy (BCG) like Q0957+561 Gl. To gauge the magnitude 
of any systematic errors thereby caused, we repeat our procedure with 42 velocity profiles from a 
set of 12 BCGs from Fisher et al. (1995), who find that the radial slopes of a p are similar to 
those of a sample of normal elliptical galaxies. Since their data do not include measurements of 
/13 and /14, we use the approximations mentioned above for v Tms and a p . As shown in Figure 
lc, the resulting BCG a p and v rms profiles are consistent to ~ 5-10% accuracy with those from 
the Bender et al. (1994) data set, although small-number statistics make the significance of the 
differences difficult to interpret. We also examine the Bender et al. (1994) data for any systematic 
correlations with other galaxy properties (the absolute magnitude Mb, the dimensionless rotation 
[£m/(6p)]*) an d the stellar projected axis ratio q*). While we find indications that "BCG-like" 
galaxies have systematically lower /14 and higher a p at large radii, they show no clear difference for 
v Tms . Although there are insufficient data available to construct BCG mean profiles, we can still 
use them to estimate the systematic corrections to our modeling results for Q0957+561 Gl (see 
§4). 
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FIG. 1. — Mean profiles computed from elliptical galaxy data. Error bars show the dispersion about the 
mean, (a) Fourth-order Gauss-Hermite moment. A solid line connects the Bender et al. (1994) data, (b) 
Velocity profiles, from the Bender et al. (1994) data, normalized to the total velocity dispersion inside the 
radius m = 0.14. Shown are both the rms velocity (solid line) and the velocity dispersion (dotted line), (c) 
Velocity profiles, from the Bender et al. (1994) (solid lines) and the BCG (dotted lines) data. For clarity, 
the scatter is not shown. 
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3. METHODS 

Schwarzschild (1979) described, and Richstone & Tremaine (1984) extended, a completely 
general method of dynamical modeling, where a galaxy is built from a library of representative 
orbits, each weighted with an occupancy number. The weights are adjusted so that the model fits 
a set of observational constraints — typically the surface brightness and line-of-sight velocities of a 
galaxy By construction, the method arrives at a solution that is a physical system of non-negative 
orbits (thereby avoiding the problems with using the Jeans equations). Unlike other common 
modeling methods, the method requires explicit knowledge of neither the integrals of motion nor 
the form of the distribution function. Recent efforts at galaxy modeling have employed sophisticated 
variants of the method that include fits to higher-order velocity moments (e.g., Rix et al. 1997; 
van der Marel et al. 1998). 

As discussed in §2.1, we adopt the SPLS family of density profiles for our models, along with 
the Hernquist (1990) mass model. The initial radii rofc of the orbits are logarithmically spaced in 
ro, and the energy of each orbit is selected to correspond to that of a circular orbit at this 
radius, <J?(ro) + v 2 (ro)/2. For a singular isothermal potential, the spacing is uniform in energy. 
The angular momentum L/, of the orbit is selected randomly from the range [0,L max ], where 
L mSLX = rov c (ro). This procedure ensures dense, uniform coverage of the (E,L) phase space. The 
model observables y m are given by the orbit weights w and the orbit projection "kernels" K m , which 
are averaged over time and over all spherical-polar viewing angles (0,<p): yf 1 = J2k w k(Kik)t,9,<j>- 
For example, the kernel for the angle-averaged surface density of an orbit at radius r is given by 
(KW^oj = (27rrvV 2 - R 2 )' 1 . The orbit is then run forward in time for one radial period T r , and 
the final kernel is found by averaging over time: {K) t fi,4> = T" 1 J Tr (K)g j( p(t)dt. To calculate the 
Gauss-Hermite velocity moments, we calculate the LOSVD £(v p ) in 41 velocity bins from v p = 
to f max , where the maximum velocity t> max is given by the largest velocity attainable on a radial 
orbit by the highest-energy orbit. We then perform a nonlinear least-squares fit to find 7 and a p , 
and then use equation (1) to find h\. 

Since the model is typically underconstrained, we fit the model observables y m to the data 
y d using the statistic \ 2 = J^iiuf 1 ~ yt) 2 / a ii while optimizing a smoothing function, the entropy 
S = J2k w k^ nw k- This corresponds to minimizing the function / = \ 2 + where the Lagrangian 
multiplier A is a smoothing factor. During our modeling runs, we reduce A slowly from 1 to 10~ 5 
to arrive at the limiting case where no smoothing is imposed. Since the velocities scale linearly 
with the total mass of the galaxy, we can leave the mass dispersion parameter (oq = l-nGp^r 2 for 
the SPLS models) free to vary in the fit. We can also enforce isotropy by minimizing the function 
fi = J2i[( uv r)i ~ ( uv t )i\ 2 /[{ vv2 )i + { uv t)i\ 2 ■ We use the conjugate gradient method (Press et al. 
1992), with first and second derivative information, to minimize /. 

We test our methods on a self-consistent (constant mass-to-light ratio) isotropic Hernquist 
(1990) galaxy, with mass density profile p(r) = Moa(27rr) _1 (r + a) -3 . We fit the exact analytic 
stellar surface brightness profile I(R) and the projected stellar velocity dispersion profile a p (R). 
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Both profiles are measured in 21 annuli from R = to R = 16a, and are assigned 10% measurement 
errors. There are 2000 orbits, spaced with initial radii from 0.07a to 221a (R e g ~ 1.8a), resulting 
in a radial coverage from to 442a. There are a variety of solutions consistent with the data, 
including the self-consistent isotropic solution. Although the more anisotropic solutions generally 
show more pronounced deviations from Gaussianity in their LOSVDs, isotropy does not necessarily 
imply hi = 0, and vice- versa (see Figure 2). 



E/|* | 
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FIG. 2. — Isotropic model solution to simulated Hernquist (1990) galaxy data. Upper left: Model (dotted 
lines) and data (solid lines), for surface brightness (with error bars) and luminosity density. Lower left: 
Same, for projected stellar velocity dispersion. Upper right: Orbit weights, in energy-angular momentum 
phase space. The area of each square is proportional to the logarithm of the orbit's weight. The bottom axis 
shows the orbit's initial radius, and the top axis its energy. Lower right: Anisotropy parameter (solid line) 
and fourth-order Gauss-Hermite velocity moment (dotted line), as a function of radius. The break radius is 
a = 1 (R cS ~ 1.8). 
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4. Q0957+561 RESULTS 

We next model the galaxy Gl in Q0957+561, fitting only the measured data: I(R) and central 
<7 P (see §2.1). We use the best-fit SPLS mass model from GN, with (a = 1.92, r c = 0'.'058). There 
are 2000 orbits with initial radii from 0'.'07 to 221", resulting in radial coverage from 0" to 361". 
A wide range of solutions fits the data exactly, and we find 1-a limits on the mass dispersion 
parameter of a$ = 295^^1 km s _1 , defined by the A% 2 = 1 boundary about the minimum \ 2 ( see 
Figure 3). Such a large range of possible solutions corresponds to a 73% uncertainty in the mass of 
Gl, but the extreme solutions show radical departures from a constant velocity dispersion profile, 
and from nearly-Gaussian LOSVDs (see Figure 4). For example, a very massive solution has nearly 
circular orbits at large radii, so that there are few plunging radial orbits to produce large velocities 
at the galactic center (see Figure 6). This behavior shows up as a velocity dispersion profile that 
rises with radius, and has a negative /14 moment (flat-topped LOSVD) at large radii (see Figure 5). 

To reject such solutions, we impose our "mean profile" constraints on a p (R) and h±(R) (see 
§2.2; note that the spherical symmetry in our model implies v rms ~ <j p ). We find that the range 
of viable solutions is dramatically reduced (Fig. 3), with new 1-a limits on <to of 280l£jj km s -1 , 
and 2-cr limits of 280±34 km s 1 . Some of these solutions may not be dynamically stable, but 
incorporating stability criteria into our model fitting is beyond the scope of this project. 




ct (km s ') 

FIG. 3. — Chi-square of the best-fit lens mass model to the Q0957+561 Gl data (surface brightness and 
central projected velocity dispersion), as a function of the mass dispersion (solid line). Also shown are the 
fits including the mean <J P (R) and hi(R) constraints (dotted line). 



-10- 



Since the reduced range of the solutions depends on the application of our mean profile con- 
straints, we next investigate the sensitivity of the results to our main concerns about these con- 
straints. First, we remove the constraints on /14 at large radii (R > 1.8i? c ff), where the reliability 
of the local galaxy data is questionable. This increases the best-fit a$ and its uncertainty by less 
than 1 km s^ 1 . Second, we estimate the effect of a systematic bias. Since we found in §2.2 that 
BCGs may have systematically lower /14 at large radii, we examine an extreme case where we set 
h A = -0.06 for R > V!2. We find new 1-a limits of a = 287l{| km s _1 , indicating that the 
systematic correction for galaxy type would increase <Tq (and Hq) by at most 5%. In summary, we 
find l-o" limits on 00 of 280^26 ^ m s_1 ' w hich corresponds to a 16% uncertainty in the mass of Gl. 
Part of the uncertainty in gq is due to the measurement error of the central a p (12 km s" 1 ), while 
part is due to the systematic uncertainty in converting from a p to 00 (~ 19 km s _1 ). This contrasts 
strongly with the conversion from GN, which implies (by assuming the anisotropy j3(r) to be con- 
stant and near zero) 1-a limits of 00 = 290^{g km s _1 . Their reported systematic uncertainty of 2 
km s _1 is clearly underestimated, given that our best dynamical model is systematically different 
from theirs by a factor five times larger. 
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FIG. 4. — Solutions for Q0957+561 Gl, for several mass dispersions. The solutions are all acceptable 
(A% 2 < 1), given only the surface brightness and central velocity dispersion data. Upper left: Surface 
brightness profile, where the error bars show data from Bernstein et al. (1993, 1997). Lower left: Fourth- 
order Gauss-Hcrmitc moment profile, with the "mean profile" shown by error bars (see §2.2). Upper right: 
Velocity dispersion profile, with the mean profile shown by error bars. Lower right: Anisotropy profile. 
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FIG. 5. — Line-of-sight velocity distributions, for the same solutions, at several radii. 
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FIG. 6. — Orbit weights for the same solutions, plotted in energy- angular momentum phase space. The 
area of each square is proportional to the logarithm of the orbit's weight. The top axis shows the initial 
radius. 
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We next compare our results with the Q0957+561 lens models of GN to produce new bounds 
on the Hubble constant Hq. Given their best fit mass model (a = 1.92, r c = 0'.'058), GN found l-a 
limits on the deflection parameter «e of 2"40 ± 0"07. As mentioned in §1, the conversion of «e to 
the Gl mass dispersion do is subject to a well-known degeneracy between the galaxy mass and the 
cluster mass: one can add a cluster with a convergence k, and decrease the galaxy mass by the 
factor (1 — k). Since this degeneracy affects none of the image observables but the time delay, Hq 
is systematically uncertain by the same factor (1 — k). Thus, with the time delay measurement of 
Kundic et al. (1997), and density parameter Qq = 1, the GN results imply do = (324 ± 4)^1 — k 
km s _1 and Hq = (82 ± 2)(1 — k) km s _1 Mpc -1 . Our dynamical constraints on o"o therefore 
put l-a constraints on k of 0.25to'io> implying Hq = 6ll^ km s _1 Mpc -1 ; and 2-a constraints 
of k = 0.25t°;{g, implying Hq = 61+^ km s" 1 Mpc" 1 (see Figure 7). Note that Tonry & Franx 
(1998) measured a central velocity dispersion for Gl of a p = 288 ± 9 km s -1 ; although their 
aperture is different from that used in our dynamical models, their measurement would roughly 
imply k ~ 0.20 ± 0.12 and Hq ~ 65 ± 10 km s" 1 Mpc" 1 . 

For comparison, we examine complementary studies of this system. Fischer et al. (1997) used 
the weak lensing of background galaxies to determine the surface density X(-R) of the cluster. Using 
their parameterized model fit to £ and their stated uncertainties, we find l-a constraints on the 
convergence k of 0.19lg q|, and 2-a bounds of (0.06-1.00). Similarly, Kundic et al. (1997) estimated 
k = 0.22 ± 0.14 (2 a) from this data, but they neglected the uncertainties in the cluster position. 
Note that these estimates of k were derived with the cluster centered on the galaxy Gl, rather 
than at the real mass centroid of the cluster, and also that the proximity of the cluster to Gl 
invalidates the description of the potential as simply a convergence k and a shear 7 (see Kochanek 
1991). Kundic et al. (1997) also warn that an error in the assumed mean redshift of the background 
galaxies can affect the derived Hq significantly. With these caveats in mind, we find that the l-a 
values for k from the Fischer et al. (1997) models imply 00 = 291^59 km s _1 for the best-fit GN 
model, and Hq = 66^27 km s _1 Mpc -1 (see Fig. 7). In a complementary study, Chartas et al. 
(1998) determined the mass of the cluster from its gaseous X-ray emission. Their results imply 
k = 0.11 ± 0.04, a = 306 ± 17 km s -1 , and Hq = 73 ± 6 km s" 1 Mpc" 1 (2 a), but they include no 
estimate of their systematic uncertainties. 

Both of these direct constraints on k are consistent with the constraints implied by our dynam- 
ical models of Gl, and could in principle be combined with them to further constrain Hq, though 
they are currently too uncertain to add any useful information. So we find that the GN lens model, 
combined with the dispersion measurement of Gl, permits a determination of Hq to 16% accuracy. 
Although better lens models may significantly shift the best-fit Hq value, they will have little effect 
on the uncertainties, which are dominated by the uncertainty in k; we have already assumed a 
perfect determination of the parameters (a, r c ) from the GN models, and a perfect determination 
of ctE would tighten the l-a limits on Hq by only ~ 0.1 km s _1 Mpc -1 . To find Hq to significantly 
better accuracy from this system, we would need much better direct constraints on the cluster mass 
distribution, and/or highly accurate velocity data from the galaxy at larger radii. 
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FIG. 7. — Constraints on galaxy mass dispersion ao, cluster convergence k, and Hubble constant Hq. The 
left panel shows a diagram of the independent constraints. The solid lines show 1-a bounds on ao and n, 
and the dotted lines show 2-a bounds. The curved lines show the constraints on <jq and k from the GN lens 
models. The horizontal lines show the constraints on n from Fischer et al. (1997). The vertical lines show 
our constraints on a . The approximate region of parameter space permitted at the 1-a level is indicated by 
dark shading; the 2-a region by light shading. The right panel shows a contour plot of the permitted region 
(including the GN and ctq constraints), with 1-a and 2-a bounds. The point marks the best-fit solution. 



5. PG1115+080 RESULTS 

We next model the galaxy G in PG 1115+080 to determine if the different lens mass profiles 
explored by Impey et al. (1998) are consistent with the measurements by Impey et al. (1998) of 
I(R) and by Tonry (1998) of the central a v (see §2.1). If any mass model can be ruled out, so can 
its corresponding value of Hq. As in the case of Q0957+561, a single central a p measurement will 
give us little information about the galaxy mass profile, so we again impose mean profile constraints 
on v Tms (R) and h±(R), normalized to this galaxy (see §2.2). 

We fit a singular isothermal mass model to these data, and although we have not exhaustively 
explored the entire range of possible solutions, we estimate the projected mass inside R = l'.T5 
to be M = 1.7lo!7 x lO n /i _1 M , which is consistent with the lens model's implied M = (1.25 ± 
0.02) x 10 11 h^ 1 Mq (for the case where the nearby galaxy group is modeled as a singular isothermal 
sphere). Similarly, we fit the modified Hubble model and find M = (1.6 ±0.5) x 10 11 /i~ 1 M Q , which 
is consistent with the lens model fit mass. To test the constant mass-todight ratio hypothesis, we 
fit a Hernquist (1990) model with a break radius of a = 0'.'325 (to match the measured effective 
radius of i? e ff = 0"59), and find M = l.St.Q^ x W n h~ 1 M Q , which is also consistent with the lens 
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model fit mass. Our dynamical models are also consistent with the lens model results when the 
group is modeled as a point mass, in which case M ~ 1.4 x lO n h~ 1 M Q . To compare the relative 
likelihood of the three mass models, we impose the lens model mass normalization on each of them, 
and find that our dynamical solutions are statistically indistinguishable (A% 2 < 1). 

With the measured time delay (Schechter et al. 1997 ; Barkana 1997b), for £Iq = 1 and a 
singular isothermal group model, the singular isothermal galaxy model gives Hq = 44 ± 4 km s -1 
Mpc -1 ; the modified Hubble profile model, 61 ± 5 km s -1 Mpc -1 ; and the constant M/L model, 
65 ± 5 km s -1 Mpc -1 (Impey et al. 1998). Treating the group as a point mass increases Hq by 
~ 10%. Since our stellar dynamical models do not rule out any of the lens models, a large range of 
values for Hq is still permitted by the system. As in the case of Q0957+561, stronger constraints on 
the mass distribution of the system (e.g., from LOSVD measurements of the lens galaxy, or from 
more detailed observations of the Einstein ring) will be necessary to break the model degeneracies. 



6. CONCLUSIONS 

Using very general orbit modeling methods, we have examined the uncertainty in the mass of 
the lens galaxies Q0957+561 Gl and PG 1115+080 G, given observations of their central projected 
stellar velocity dispersions a p (Falco et al. 1997; Tonry 1998). As many past studies have shown, 
such a measurement alone is inadequate to strongly constrain the galaxy's mass. In order to put 
additional realistic a priori constraints on the galaxy's properties, we have derived "mean profiles" 
of the rms projected velocity v rms (R) and the fourth-order Gauss-Hermite moment h^(R) from 
a large sample of nearby elliptical galaxies. These mean profiles prove to be remarkably self- 
similar, even over a large range of galaxy types. This universality is not too surprising, given 
the homology of early-type galaxies implied by the existence of the fundamental plane — an even 
stronger implication if the central kinetic energy is considered instead of the velocity dispersion 
(e.g., Busarello et al. 1997). 

For Q0957+561 Gl, given the best-fit SPLS mass model from Grogin & Narayan (1996), with 
only the surface brightness profile I(R) and central o~ p as constraints, we find 1-cr limits on the mass 
dispersion gq of 295^21 km s_1 - The addition of the mean profile constraints reduces this permitted 
range to 280^^6 km s_1 - I n conjunction with the GN lens model constraints, this implies a cluster 
convergence of k = 0.25lg;^, which is consistent with the constraints on k from other independent 
studies. Using the time delay measurement of Kundic et al. (1997), we find l-a limits on Hq of 
61^1! km s _1 Mpc -1 , and 2-a limits of 6ll]J km s _1 Mpc -1 . Thus, current measurements of the 
lens system Q0957+561 do not constrain Hq to better than 15%. To obtain useful limits on Hq, we 
will need better constraints on the cluster convergence and/or better velocity profile measurements 
for the galaxy Gl. We will also need substantial revision of the lens model — the GN solutions 
fit the lens data poorly, inaccurately assume a spherical galaxy, and use an oversimplified Taylor 
expansion model of the cluster. Presumably a more accurate lens model will eventually provide 
a good fit to the data and a different value for Hq (see Barkana et al. 1998 for some improved 
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models), in which case our dynamical model will still be illustrative of the systematic uncertainties 
expected. 

We have also examined the lens galaxy PG 1115+080 G, which has a total mass that is relatively 
well-determined by the lensing constraints, but a radial mass distribution that is unconstrained 
(Impey et al. 1998). As with Q0957+561 Gl, we model this galaxy by including constraints on 
I(R) and the central a p along with the mean profile constraints. We find that these constraints are 
not sufficient to break the degeneracy between the different lens models, so that a large range of 
values is still permitted for Hq (44-68 km s" 1 Mpc -1 ). Further constraints on the mass distribution 
in this system are needed. 

The stellar dynamics of gravitational lens systems show considerable promise for determining 
Hq, even if we have only central dispersions. Larger samples of dispersion measurements would be 
particularly valuable, since each lens system will have a different set of systematic uncertainties (it is 
encouraging to note that the results for Hq from Q0957+561 and PG 1115+080 are so far consistent 
with each other) . Of particular value would be results from isolated systems with lens geometries 
that probe the radial mass distribution, such as MG 1654+1346 (see Ellithorpe, Kochanek, & 
Hewitt 1996) and MG 1549+3047 (Lehar et al. 1993, 1996). The results from a large set of these 
systems could be combined using Bayesian methods to converge on a robust value for Hq (see Press 
1997), avoiding the correlated systematic uncertainties that plague "distance ladder" approaches. 
Furthermore, the independent measurement of galaxy properties from gravitational lensing and 
from stellar dynamics has an enormous potential for shedding light on the detailed structure of 
galaxies. 
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